Mon. Not. R. Astron. Soc. 000, 



(0000) 



Printed 1 February 2008 (MN IMfeX style file vl.4) 



Radio sources near the core of globular cluster 47 Tucanae 

D. McConnell 1 and J.G.Abies 2 

1 ATNF, CSIRO, Locked Bag 194, Narrabri NSW 2390, Australia 

2 Telecommunications and Industrial Physics, CSIRO, PO Box 76, Epping, NSW 2121, Australia 



o 

Q 

(N 



> 

in 

m 

(N 

on : 

ON ■ 

^ ' 

Or 

i ' 
o ■ 

a: 



Version: 1 February 2008 



ABSTRACT 

We present ATCA radio images of the globular cluster 47 Tucanae made at 1.4 and 
1.7 GHz and provide an analysis of the radio sources detected within 5' of the cluster 
centre. 11 sources are detected, most of which are clustered about the core of 47 Tuc. 
Both of the pulsars in 47 Tuc whose positions are known can be identified with sources 
in the 1.4 GHz image. The source distribution has a characteristic radius of ~100", 
larger than the 23" radius of the cluster core. We compare source positions with the 
positions of nine X-ray sources and find no correspondence. 
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1 INTRODUCTION 

The globular cluster 47 Tucanae is known to contain many 
s hort period (millisecond) radio pulsars. Robinson et al. 
(1995) reported 11 pulsars observed with the Parkes radiote- 
lescope at 640 MHz and 436 MHz. Millisecond pulsars as 
a class are intrinsically faint. Coupled with the relatively 
great distance to 47 Tuc (4.6 kpc), this means that most 
of the pulsars are difficult to detect. Robinson et al. report 
that many detections rely on the occasional flux enhance- 
ments produced by the focussing effects of diffraction and 
refraction in the irregular interstellar medium. Many of the 
pulsars show period variations characteristic of motion in 
a binary orbit. The use of pulse timing to determine pul- 
sar positions in the cluster has been made difficult by both 
their intermittent detectability and the extra complexity of 
timing analyses arising from their orbital motions. In 1995, 
Robinson et al. were able to report positions for only two of 
the pulsars — 47 Tuc C and 47 Tuc D. 

The positions of the pulsars relative to the cluster cen- 
tre are interesting because of the relation to the dynamics 
and evolutionary history of the cluster core. For example, 
pulsars C and D have a negative period derivative which is 
attributed not to an intrinsic increase in rotation rate, but to 
their accelerated motion in the cluster's gravitational field. 

Here we report radio images of 47 Tuc made at 1.4 GHz 
and 1.7 GHz with the Australia Telescope Compact Array 
(ATCA) and an attempt to seek point radio sources which 
may be identified with pulsars in the cluster. At 1.4 GHz, 
the sky away from the Galactic plane is dominated by ex- 
tragalactic sources with small angular sizes. The number 
density of these sources increases, naturally, with decreas- 
ing flux density and at flux levels expected of the radio pul- 
sars in 47 Tuc, tens of sources are expected in each ATCA 



primary beam. Several reports of studies of 1.4 GHz source 
densities are available as a guide to the expected number of 
backg round sources. In pa rticular, the "Phoenix Deep Sur- 
vey" (Hopkins et al. 1998) was conducted at 1.4 GHz with 



the ATCA and provides a useful reference. In the analysis of 
these observations of 47 Tuc we study three source charac- 
teristics which might distinguish the radio pulsars from the 
background sources. (1) Spectral index: the spectral index 
a (where S ~ u a ) of pulsar emission tends to be less than 
(a < —1.2) the typical value for extragalactic sources for 
which a ~ —0.8. (2) Flux variability: the very small angu- 
lar size of radio pulsars leads to large flux variations due to 
scintillation in the interstellar medium. (3) Spatial distribu- 
tion: background sources are distributed uniformly whereas 
the pulsars are expected near the cluster core. 

Radio images of other globu lar clusters have been re- 
ported by Kulkarni et al. (1990). They observed five clus- 
ters with the VLA and detected the pulsars in clusters M4 
and M28, one of now eight known pulsars in M15, a single 
unidentified radio source in M3 and no radio sources in M92. 



2 OBSERVATIONS AND DATA REDUCTION 

The Australia Telescope Compact Array (ATCA) was used 
to observe 47 Tuc at 1.408 GHz and 1.708 GHz. At each fre- 
quency a 128 MHz bandwidth was used, each band split into 
64 spectral channels. The cluster was observed on 1992 Jan 
24 in the 6D array configuration and on 1992 Apr 23,24,25 
in the 6C configuration. Four fields close to the cluster cen- 
tre were observed. Table 1 gives the field centres and their 
times of observation. 

The central field was observed in two separate 12-hour 
integrations which are listed individually as I and II in the 
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Table 1. Summary of observations. Six images were made at each frequency, five from each 12-hour integration and a sixth, the mosaic, 
as a combination of I-V. The sensitivity of each observation and the mosaic is indicated by the root mean square value of pixel brightness 
(ex), expressed as /ijy/beam, and measured over the central 10' of the cluster. 



Observation a S Date Time 01.4 Cl.7 

(J2000) (J2000) (UT) (UT) (/LtJy) (//Jy) 
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table. Fields III IV and V were offset from the cluster centre 
by 15', less than the ~17' to the half-sensitivity point of the 
ATCA primary beam at 1.4 GHz. Thus all five integrations 
individually had good sensitivity to emission from the clus- 
ter centre, and in combination provided a sensitive image 
over a 35' region centred on the cluster. 

The miriad data reduction package was used to cali- 
brate the data and form images. The antenna gains were 
calibrated using brief observations of the source B2353-686 
(5*1.4 = 1.05 Jy, 5*1.7 = 1.00 Jy) and the flux scale was 
referred to the ATCA primary flux calibrator B1934-638 
(5i.4 = 14.9 Jy; Si. 7 = 14.0 Jy). Images were formed for 
each integration at both frequencies with reference posi- 
tion (tangent point) at a(2000) = 00:24: 05° 83, (5(2000) = 
-72:04:5l'.'4, the nominal cluster centre (Guhathakurta et 
al. 199^). Visibility data from each 2 MHz spectral channel 
were gridded separately allowing the full primary beam to 
be imaged with very little bandwidth smearing. An imaging 
cell size of l'.'5 was used to provide adequate sampling of the 
synthesized beam. The image for each integration was made 
2041 pixels square (51'). The final image at each observing 
frequency was formed as a linear combination of the images 
from each field. In the remainder of this paper we discuss a 
central region of the image of diameter 10' and centred on 
the cluster. 



3 IMAGE ANALYSIS AND RESULTS 
3.1 The image noise statistics 

Five sub-images, corresponding to the integrations listed in 
Table 1, were made at each frequency. A sixth, the final (mo- 
saic) image, was formed by linear combination of the five. To 
allow direct comparison, the five sub-images were corrected 
for the sensitivity variation across the primary beam. Table 
1 gives the noise in each sub-image in the central 10' of the 
cluster as the root-mean-square of their pixel values (<r). The 
higher noise levels in sub-images IV and V are possibly due 
the effects of imperfect antenna pointing during the integra- 
tion and the presence of relatively strong sources (92 mjy 
and 41 mjy) near the halfpower levels of those fields. The 
size (FWHM) of the synthesized beam in the combined im- 
ages was 5'.'7 x 5'.'2 at 1.4 GHz and 4'.'5 x 4'.'0 at 1.7 GHz. 



3.2 Radio sources 

Examination of the images shows a number of point sources. 
Sources within a 5' radius of the cluster centre were anal- 
ysed. All sources with a peak brightness exceeding 4.5<j in 
the 1.4 GHz mosaic were selected. Their positions and fluxes 
were estimated by fitting each with a gaussian function 
with major and minor axes equal to those of the synthe- 
sized beam. Because of the known flux variability of pul- 
sars in the cluster, each of the sub-images (I - V) were 
searched for sources coincident in two sub-images but that 
were not strong or persistent enough to rise above the detec- 
tion threshhold in the mosaic. No such sources were found. 
All sources detected at 1.7 GHz were also visible at 1.4 GHz. 
Table 2 lists the sources detected in the 1.4 GHz image. The 
1.7 GHz image was analysed in the same way and 1.7 GHz 
fluxes are reported for each source in Table 2. Two sources 
(2, 4) had a 1.7 GHz brightness below the 4.5a", threshhold, 
but for each a successful fit to a model gaussian source was 
made, yielding 1.7 GHz fluxes of 165/xJyand 177^tJy(3.6<7 
and 3.8a) respectively. Two further sources (9, 11) had no 
detectable 1.7 GHz emission and upper limits correspond- 
ing to 3cr are quoted in Table 2. The radial distance of 
each source from the cluster centre (a(2000) = 00:24:05183, 
5(2000) = -72:04:51'.'4, Guhathakurta et al., 19920) is listed 
in the table. Table 2 also lists, in column 8, the spectral 
index of each source. Fig. 1 shows the location and relative 
flux density of the 1.4 GHz sources superposed on an optical 
image (from the Digital Sky Survey) of the cluster. 

The uncertainties in source positions derived from the 
numerical model fitting do not account for all contributing 
effects and are typically too small by a factor of ten. Two 
alternate methods were used to estimate the uncertainties in 
source position. Several sources are detected in several of the 
five sub-images which allowed the dispersion of measured po- 
sitions to be determined directly as their standard deviation 
A9m- This quantity was compared with A9 C ~ FWHM/snr 
where FWHM ~ 5" is the synthesized beam size and the 
signal-to-noise ratio snr — S1.4/& — 5i.4/42/xJy. For the 
sources with a number of independent position measure- 
ments, A8 m and A9 C , differ by less than 20% and so A6 C 
was used for all sources. The uncertainties in right ascen- 



* Note that an inde pendent determination of the cluster centre 
(Calzetti et al. 1998) differs by 1'.'4 which is not significant in this 



analysis. 
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Figure 1. The positions of sources detected in the 1.4 GHz image of 47 Tucanae shown against an optical image (from the Digital 
Sky Survey) of the cluster. The large circle marks the boundary of the region analysed. Each point source is marked as a filled circle 
with size proportional to the logarithm of the source flux Si. 4. The faintest sources indicated have Si. 4 > 190/xJy. Two extended sources 
of emission are indicated near a(2000) = 00:24:47, 5(2000) = -71:59:15 and a(2000) = 00:25:04, <5(2000) = -72:02:36 by a contour at 
220/^Jy/beam. These regions have integrated 1.4 GHz fluxes of approximately 60 mjy and 6.5 mjy respectively. The small central circle 
marks the position and size (r c ~ 23") of the cluster core. The locations of pulsars C and D as determined from pulse timing are indicated 
by crosses (pulsar C to the west). 



sion and declination are assumed independent and quoted 
in Table 2 as A a — V2A6 C and AS = ^2A0 C . 



3.3 Source distribution 

From Fig. 1 it is evident that the radio sources are clustered 
near the image centre. Our analysis of source positions is 
illustrated in Fig. 2 which shows the cumulative source count 
as a function of distance from the cluster centre. To quantify 
the size of the central group, the radial density variation 
N(r) has been modelled by a gaussian function added to a 
uniformly distributed background distribution 

N{r) =ae- {r/b>2 +c 

The fit given by a = 0.96 arcmin -2 , b = 91", c = 
0.035 arcmin" 2 is plotted on Fig. 2. The value of the param- 
eter c (which corresponds to 2.7 sources in the analysed area) 



can be compared with the expected density of backgr ound 
sources. From an analysis of the Phoenix Deep Survey ( [Hop- 



kins et al. 1998), taking account for extended sources falling 



below the brightness detection limit and assuming a Pois- 
son distribution of source counts in a given area of sky, 
we find an expected value for c of 0.10 ± 0.04 arcmin" 2 or 
7.8 ± 2.8 sources over the 10' field. The number of sources 
detected in the field is similar to the expected number of 
background sources. Note however that 7 of the observed 
sources lie within 2' of the cluster centre, compared with 
the 1.3 ± 1.1 sources expected from the Phoenix Deep Sur- 
vey source counts. 



3.4 Source variability 

Obse rvations of the pulsed emission from the pulsars in 47 
Tuc ( [Robinson et al. 1995|) show their flux density to vary 
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Table 2. Point sources within 5' of the cluster centre. Sources with positions corresponding to known timing positions of pulsars are 
indicated by the pulsar letter designation in column 2. Column 5 gives the distance of each source from the cluster centre o(J2000) = 
00:24:05?83, 5(J2000) = -72:04:51"4. The uncertainties given for source positions (columns 3, 4) relate to the last digit or digits quoted 
and were determined by the technique described in the text. 
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Figure 2. Cumulative count of source number inside radius 
r. The curved line indicates the cumulative source count from 
model density distribution N(r) = ae~( r < /!> ) 4- c with a = 
0.96 arcmin -2 , b = 91", c = 0.035 arcmin -2 . 



strongly, presumably due to interstellar scintillation. How- 
ever, most extragalactic radio sources have a larger angular 
extent (while still appearing point-like with the ATCA) and 
so do not scintillate, or do so only weakly. The observations 
described here allow us to make several flux measurements 
of some of the 47 Tuc sources. Each sub-image was exam- 
ined at every source location listed in Table 2. A source was 
deemed visible in a sub-image if its flux exceeded 3 times 
the image rms from Table 1. Seven of the 11 sources had 
measureable flux by this criterion in three or more of the 
sub-images. Their fluxes are plotted in Fig. 3. 



4 DISCUSSION 

We have presented in Table 2 a list of 11 radio sources within 
5' of the centre of 47 Tuc. The positions of two of these, 
4 and 8, ca n be identified with pulsars 47 Tuc C and D 
respectively (Robinson et al. 1995). These three sources have 
steep spectra, 04 = —2.9 and as = —1.4 and varying flux 
denstities as indicated in Fig. 3. Robinson et. al. report the 
spectral indices of 47 Tuc C (0.0) and D (-1.4) measured 
between 436 MHz and 640 MHz, suggesting that 47 Tuc D 
has a straight spectrum with slope -1.4 from 436 MHz to 
1.7 GHz but that the spectrum of 47 Tuc C steepens above 
640 MHz. 

Source numbers 1 and 10 in Table 2 are likely to be 
background sources — both are distant from the cluster cen- 
tre and have spectral indices typical of extragalactic objects. 
Source 1 shows little flux variation between the five obser- 
vations. Source 3 also has a flat spectrum and is probably 
extragalactic, although some flux variability is suggested. 
Sources 2 and 5 both exhibit flux variations, 5 strongly so, 
and have very steep spectra and so it seems likely that both 
are pulsars. The great distance of source 2 from the clus- 
ter centre (248") makes it an interesting pulsar candidate. 
Source 5 has a steep spectrum and variability characteristic 
of a pulsar and is stronger than sources 4 and 8 (pulsars 
C and D), but the pulse timing analyses of Robinson et 
al. (1995) did not yield pulsar parameters for its location. 
Perhaps source 5 is a pulsar with properties that make the 
timing analysis difficult or that cause its pulses to be un- 
detectable, such as a pulse period outside the range of the 
original observations, an extraordinary binary motion or a 
dense local environment strongly scattering the pulsed radi- 
ation. The case for the apparently constant source 6 being 
a pulsar is not strong, although it has a spectrum a little 
steeper than expected for background objects. Sources 7, 
9 and 11 are faint, making variability and spectral shape 
difficult to measure. We note a similarity in the pattern of 
intensity variations between sources 4 and 8 (Fig. 3). We 
can find no instrumental effect which could lead to such ap- 
parent variations that would not similarly effect the other 
sources in the field. On the other hand, variability induced 
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Figure 3. 1.4 GHz fluxes of a selection of sources (identified with 
entries in Table 2 by the number in the upper right) in the field. 
Up to five flux values are plotted for each source, corresponding 
to the five independent observations listed in Table 1. The hori- 
zontal broken line indicates the source flux measured in the final 
mosaiced image. The error bars on each point indicate the rms of 
pixel brightness in each case. 



by scintillation in sources as separated as these is difficult 
to understand. 

Verbunt and Hasinger (1998) report the positions of a 
number of X-ray sources in 47 Tuc, nine of which lie close to 
the cluster core. Fig. 4 shows their locations relative to the 
radio source positions and the locations of pulsars C and D. 
It is evident that the X-ray sources belong to two popula- 
tions, those within the cluster core, and those represented 
by X4, X6, Xll and X13 which have a radial distribution 




-72 06 30 - 

h 24 m 30 s 20 s 10 s 24 m B 50 s 40 s 

RA (J2000) 

Figure 4. The inner part of the imaged field. The radio emission 
is shown as the 1.4 GHz brightness contours at 180, 350 and 700 
/ijy/beam, and the sources identified by their sequence number 
from Table 2. The size of the synthesized beam at 1.4 GHz is 
indicated at lower left. The cluster core is marked as the open 
circle in the centre. Locations of the nine X-ray sources reported 
by Verbunt and Hasinger (1998) are shown as open squares and 
labelled with numbers from Vebunt and Hasinger's Table 2. 



similar to that of the radio sources. Verbunt and Hasinger 
are unable to provide a positive identification for any of 
them, but point out that their luminosities are consistent 
with them being soft X-ray transients, cataclysmic variables 
or recycled pulsars. The first epoch of the X-ray detections 
reported by Verbunt and Hasinger is close to the date of our 
most sensitive observations (1992 April). We note that none 
of the X-ray sources coincide with any of the radio sources. 

We have presented a radio image of the central region of 
47 Tuc and shown a concentration of radio sources close to 
the cluster core. Although a direct correspondence between 
radio sources and pulsars is possible in only two cases, it 
seems likely that several of the other detected sources are 
also pulsars. Thus the distribution of sources in the image is 
likely to be representative of the spatial distribution of pul- 
sars in the cluster. We can now differentiate several classes 
of object, each with its own characteristic radial distribu- 
tion. There is a class of X-ray object represented by sources 
X5, X7, X9, X10 and X19 in the catalogue of Verbunt and 
Hasinger (1998) which are confined to the stellar core of the 
cluster. There is a second class of X-ray object (represented 
by X4, X6, Xll, X13, Verbunt and Hasinger) which lie typ- 
ically 2-5 core radii from the cluster centre. Finally there 
are the radio sources (possibly pulsars) whose distribution is 
similar to that of the outer X-ray class, but which may have 
some members very far (~10 core radii) from the cluster 
core. 
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